Heavy-quarkonium hadron cross section in QCD at leading twist 
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We compute the total cross section of a iieavy quarkonium on a liadron target in 
leading twist QCD, including target mass corrections. Our method relies on the an- 
alytical continuation of the operator product expansion of the scattering amplitude, 



' obtained long ago by Bhanot and Peskin. The cross section has a simple partonic 



form, which allows us to investigate the phenomenology of J/ip and T dissociation 



by both pions and protons. 
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(N ; I. INTRODUCTION 

(N 

^ ' It has been conjectured that in ultra-relativistic collisions between heavy ions a quark 

. gluon plasma (QGP) is formed for a very short period of time. Later it disintegrates into the 
Qs^ \ hadrons which are finally seen in the detectors. It is the challenge of the present experiments 
O \ at the CERN-SPS and at the Relativistic Heavy Ion Collider (RHIC) to find observables 
^ ■ which unambiguously signal the formation of the QGP. 

One of the most promising suggestions advanced so far is the suppression of charmonium 
production in these collisions. It has been argued by Matsui and Satz |I| that the interaction 
^ \ between a heavy quark Q and antiquark Q is screened in a QGP, and consequently the QQ 
bound pairs may not survive in this environment. The problem with this signal is that all 
Oh. other possible suppression mechanisms have to be well understood. 

r-| ! Whereas the general features of charmonium production in proton-proton collisions seem 

to be under control, already in proton- nucleus reactions the suppression is not well under- 
stood so far. Only recently, data have been published which show a different suppression of 
J/ip and ip' 0], and hence give the first hints that the J/ip is formed inside the nucleus. 
5^ \ In heavy ion collisions the situation is even more difficult. There many particles are pro- 

duced which possibly collide with a charmonium and may cause an observable suppression 
even if a QGP is not formed at all. In order to quantify such a suppression it is necessary 
to know the strength of these interactions. Experimentally the charmonium-hadron disso- 
ciation cross sections are not accessible. Therefore one has to rely on theoretical estimates. 
Three kinds of approaches have been advanced in the past. 

The first approach is based on twist expansion techniques well known from deep inelastic 
scattering studies. It has been launched by Bhanot and Peskin 0, |] and was explicitly used 
in Ref. for ip-p in the approximation of a vanishing proton mass. The whole approach 
gives a correct approximation of QCD provided that the heavy quark mass is large enough. 
It has thus the very advantage to be a well defined approximation scheme of the underlying 
theory. For realistic systems, such as charmonia and bottomonia, power corrections are 
presumably non negligible, but some of them, namely finite target mass corrections, can be 
incorporated in a systematic way. 

A second attempt is formulated within a constituent quark model. In an early study 
Martins et al. H have calculated & J /ip dissociation cross section a^-j^ by vr's of up to 7 mb 
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at ^/s = 4 GeV, i.e., 0.8 GeV above mj/^ + 771,^. This value has been reduced to about 1 mb 
for the same energy in a more recent study by Wong et al. who used parameters adjusted 
to other elementary reactions. 

A third approach is based on hadronic degrees of freedom. Invoking a local U(4) sym- 
metry and employing pseudoscalar-pseudoscalar-vector coupling, Matinyan and Miiller 
investigated the dissociation of the J/ ip by exchange of a D or Z) meson. Employing vector 
dominance to determine the coupling constants they arrive at a^-,^ ~ 0.3 mb for ^/s = 4 GeV. 
Later, Haglin |^ included four point interactions and a three vector-meson coupling and ob- 
tained a much larger cross section because the large suppression of the cross section due 
to the D-meson propagator is not present in the contact terms. Recently Lin and Ko |^ 
modified the details of this approach and included form factors. Depending on the form 
factor assumed they get cx^^ between 4 and 25 mb at y/s = 4 GeV. 

The purpose of this article is to extend the work Ref. |^ in three different directions. First 
of all we include systematically the masses of the scattering partners using a method known 
from deep inelastic scattering studies. This allows for the calculation of the dissociation cross 
sections close to threshold where it is most relevant for the question at hand. Second the 
cross section of Ref. is derived in a different, more direct fashion, again in analogy to the 
calculation of the forward Compton scattering amplitude in the operator product expansion. 
This cross section has a simple partonic expression, even when target mass corrections are 
included. We also explicitly derive how the reaction threshold is shifted by target mass 
corrections and how the cross section is modified in the vicinity of threshold. Third the 
calculation is extended towards other hadrons h and towards bottomonium which becomes 
an observable particle in the upcoming experiments at the RHIC and at the Large Hadron 
Collider at CERN. 

Target mass correction for J/ip-p cross section has been examined in the framework of 
Bhanot and Peskin in Ref. |Ty]. These authors obtained their results in the form of sum 
rules. In principle these sum rules contain all of the above mentioned aspects but none is 
made explicit. Further we were neither capable to reproduce their exact expressions for the 
sum rules nor to find the trend they mention for the correction. We will clarify in the course 
of this study where we disagree. 



II. DERIVATION 

In this section we generalize the expression obtained by Bhanot and Peskin for the 
total cross section of a heavy quarkonium $ with a target hadron h by including finite target 
mass terms. The proposed analysis is close to that performed in the context of deep inelastic 
scattering 0, O. 



A. Short review of the framework 

Let us first collect the material we need from Refs. p, ^. We want to compute the $-/t 
total cross section^. Our starting point is the expression for the forward elastic scattering 
amplitude, M.q>h- This amplitude depends on energy and it is convenient to express it in 



Throughout this paper we restrict ourselves to spin-averaged cross sections without further mention. 
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terms of 



{K + pf -M'^ 



m 



2 

2M 



where and p are the $ and h respective 4-momenta and M and rrih their respective 
masses. We note that A is the hadron energy in the $ rest frame. Via the optical theorem, 
the forward scattering amphtude leads to the ^-h total cross section 

o-<i./i(A) = ^ lmM,i,h{X)- 

Notice that we use the same definition for M. as in Ref. [Q. 

In QCD, in the limit of a large heavy quark mass, the scattering amplitude has a twist 
expansion. In the $ rest frame the Leading Twist (LT) contribution is Q 

-MiT(A) = alel E d,,e-,^\h\F'^{^D'r-'F^h). (1) 
fc>i 

ao and eo are, respectively, the Bohr radius and the Rydberg energy for the QQ system. 

The above formula displays the factorization of the process in terms of hard coefficients 
dn and soft matrix elements. Both should be evaluated at a factorization (and renormaliza- 
tion) scale /i, to be chosen to minimize the influence of neglected higher order perturbative 
corrections. It is argued in Ref. that /i ~ eo, though a precise determination fulfilling the 
latter requirement would need a complete one-loop computation. In the phenomenological 
study, we will quantify the consequences of this scale uncertainty. 

The coefficients dn correspond to matrix elements of definite operators evaluated in the 
$-state. These are computable in perturbative QCD and have been made explicit in Ref. 
for IS and 2S <l>-states to leading order in the coupling and to leading order in where 
Nc is the number of colors. For IS-state^ they read 

16^ 

rf„ = ^ 5(^ + 5/2,5/2), 

where -B(/U, v) is the Euler beta function. For later convenience we remark that dn can be 
expressed as the n-th moment of a given function / through 

1 dx 



dn / f ("^ 

Jo X 



with 

For every k in Eq. (|l|) a gluon twist-2 operator evaluated in the hadron h {spin- averaged) 
state also appears. Each of these matrix elements is a traceless fully symmetric rank 2k ten- 
sor built from the hadron momentum p^^. It turns out that a tensor having these properties 



The present consideration applies to 2S-state as well but is left out until Appendix 
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is necessarily proportional to [|TT 



X 



g®V® ■ ■ ■ ®P 



.perm. 



2k-2j 



where the tensor in the jth term on the right-hand side is the sum of the {2k)\ / [2^ j\{2k — 2j)\] 
distinct tensors one can construct by multiplying j g^^^s and {2k — 2j) p'^'s. The matrix 
element needed in Eq. (|l]) is therefore proportional to n° "'^(p) and writes 



{hip)\F'''{tDy-'F,'\h{p)) = A2,n°-°(p) 



(2) 



This set of matrix elements is related to the unpolarized gluon density G in the hadron 
target. One can see this from the matrix element definition of G, see, e.g., [|13|, which in the 
light cone gauge A'^ = reads 



xG{x) 



p^ 



dy 

271 



.e^-p'-y- {h{p)\F+^{0)F^ + {y-)\hip)). 



In the parton model the argument of G, i.e., x, is interpreted as the fraction of the hadron 
light cone momentum p'^ carried by the gluon. Taking the nth moment of G{x), we get 

l'^x-G{x) = j^{h{p)\F+''{zd^r-'F^ + \h{p)). 

Since = in the = gauge we recognize 

{h{p)\F+''{zD+)''-'F, + \h{p)) = A2kU+-+{p) 



A2k (p^ 



\2k 



I.e., 



A^ = f^—x^G{x). 
Jo X 



B. Massless target 

In the present subsection we want to illustrate the general method in the case of vanishing 
target mass, rrih = 0. Then Eq. (0) simplifies to 

and Eq. (0) thus becomes 

AliT(A) = alel Y: d2kA2k{X/eor. (3) 

k>l 
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It is useful to study the scattering amplitude A^#/i throughout the complex plane of 
the energy. To avoid confusion, we will from now on reserve the notation A to real values 
and define A as the extension of A to complex values. Using the d'Alembert criterion one 
easily checks that the convergence radius of the power series now considered with the 
complex argument A, is equal to eo- As extensively discussed in Ref. the twist expansion 
of the scattering amplitude provides an expression for A^^P(A) in the unphysical region 
of energies. Since we are interested in physical energies we have to perform an analytic 
continuation of the power series. 

(LT) 

Before doing so let us first elaborate on possible differences between A4.)^h the full 
scattering amplitude J^<s>h, i-e., including higher twist terms. It turns out from the analysis 
below that the singularities (branch points) of on the boundary of the convergence 

disk of the power series lie at A = ±eo. This is not what is expected for the locations 
of the (first) branch points of i.e., the locations of the thresholds for both reactions 

— > X and ^ + h ^ X , which are A = ±m/i. The technical reason for the difference is, of 
course, that rrih occurs nowhere in Eq. (^. In the twist expansion approach the locations of 
the singularities of the full scattering amplitude may be affected by higher twist corrections. 
We verify that this is the case for the above reactions by taking into account elastic unitarity 
corrections for which the thresholds are clearly located at A = ±m/i.^ In the next section 
we shall see how rrih 7^ corrections affect the locations of the LT reaction thresholds. 

Since the convergence radius of the power series (^ is non zero the LT amplitude"^ can 
be unambiguously determined throughout the A complex plane from the sole knowledge of 
Eq. (0), using the Mellin transform machinery. We first remark that M„ = (i„y4„ being a 
product of moments one can express 



with 



/' — x"/i(x), (4) 

JO X 



h{x) = G(g)f{x), 
and the convolution product defined as 

G®/(x)= f''^G{y)f{x/y). 
Jx y 

Now, plugging Eq. (1) in Eq. (0) for | A| < eo, freely interchanging summation with integration 
and summing the ensuing geometrical series, one finds 



ox 1 — x^(A/eo) 



2 ' 



The key point is that this integral representation can be extended throughout the entire 
complex plane except for the two branch points A = ±eo. The analytic continuation of 



^ Let us note in passing one important phenomenological consequence of the LT analysis : what we call 
the total (LT) cross section does not in fact include processes such as the elastic one. The word "total" 
is thus misleading, at least in the threshold region. 

^ From now on A4 always refers to the leading twist part of the forward elastic scattering amplitude 
and we drop the indices (LT) and $ h for simplicity. 
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Eq. to energies A > eo is then easily derived 



and its imaginary part given by 

ImA^(A) = ^[M{\ + ie) - M{\ - ie)] = ^a^el h{e^/X). 
Putting all things together, one gets 

ImA^(A) = T^oeo^/^^/xG(x)^-3^^-^. (6) 

Dividing Eq. by A (the flux factor when nih = 0) we recover the partonic expression 
of the ^-h total cross section as obtained by Bhanot and Peskin within a parton model 
approach, i.e., 

o"$/i(A) = f dxG{x) a^g{xX), (7) 

with the $-gluon cross section 







corresponding to the gluon energy in the $ rest frame. Aside from its energy dependence, 
the $-gluon cross section is driven by a^eo oc asal, as expected in QCD for the interaction 
of a small color singlet dipole of size Oq. 

In this formulation one important physical aspect is made transparent the leading 

twist analysis describes the $ dissociation by gluons into a Q and a Q 

To be energetically possible the gluon energy has therefore to be larger than the QQ Coulomb 
binding energy eo. In view of the fact that the conflnement scale is small as compared to 
eo the LT analysis then provides a description of $ dissociation into open channels, e.g., 
Qq + Qq. Let us emphasize that this dissociation is precisely the process of interest for the 
question of $ suppression in heavy ion collisions. 

An important aspect for the phenomenology of the above cross section is its limiting 
behaviors for both small and large energy regimes. These are linked to the x ^ 1 and x ^ 
behaviors of G, respectively. It is then convenient to have in mind the simple, yet standard, 
parameterization 

G{x) = A{l-xy/x^+^. (9) 

With this ansatz one can write down exact asymptotic formulas either by following the 
reasoning of Ref. |^ or by noticing that the (^-h cross section is proportional to a hypergeo- 
metric function. This is most easily done by changing variable xtot = (xA/eo — l)/(A/eo — 1) 
in Eq. (|^). Then one recognizes 

1 fi^TT 

c^(A) = -^ageoAi?(r/ + l,5/2)(A/eo-ir^/^ 



X 



{X/eoY-\F^{5 + 6, 5/2; r] + 7/2; 1 - A/eo). 
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For A in the neighborhood of eo, i.e., the rrih = threshold, the hypergeometric function 
approaches 1 and we get 

1 fi^ 77 

a<,H{X) ~ f^aleoAB{v + l,5/2) (A/eo - 1)^+^/^ (10) 



For large energies, using 



5.„ ,7., , _ B(* + 7/2,5/2),^,^ ,,5,, 



J.(^ + 6.-;, + -;l-AA„)~ ^(„/,_-/ (AA. 



one obtains 



a*.(A)~^a^eoAi?(5 + 7/2,5/2)(^A^ • (H) 

The high energy cross section is primarily geometrical (remember a^eo oc a^ag). In addition 
to this simple behavior, there is a non trivial energy dependence coming from the small x 
behavior of the gluon density. 

For phenomenological investigations we shall also use slightly more involved forms for 
G{x) as obtained in parton distribution function studies. In this case the connection to 
2F1 is lost. One may, however, derive similar asymptotic expressions by first expanding the 
gluon distribution either in the neighborhood of 1 or 0. 



C. Massive target 

Having illustrated the method for the case rrih = 0, we now turn to the general case 
nih 7^ 0. Plugging Eq. (H) into Eq. ([l|) leads to 

k>l 

i2k-j)l ,2fc_2,7 

with = c?2fc^2fc. We thus get an amplitude which may be considered as a double power 
series in A and rrih. The study of this double series with complex arguments A — A and 
m|/(4eQ) — > z shows that it is absolutely convergent for |A/eo| + \z\ < 1. In this domain, 
defining k' = k — j we may rewrite the series as 



{2k' + j)\ 



X 



The second term on the right-hand side corresponds to the power-series expansion of the 
scattering amplitude A4 of Section II B evaluated at the complex plane location A = imh/2: 
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From the representation (j^) of the scattering amphtude we immediately see that M.{imh/2) 
is well defined and real for every (real) irih. We thus ignore this term in the following since 
it does not contribute to the total cross section at leading twist. 

Let us now concentrate on the first term on the right-hand side of Eq. (p!2D. We use 
the same reasoning as in Section LIB considering now the double series (12) with complex 



arguments A and z instead of A and m^/(4eQ). Expressing first M„ as the n-th moment of 
h = G ® f one can write 



{2k' + j)\ , , . /■! dx 2fc', , . 



j\{2k') 



-x' zy 



In the convergence domain of the double series, the series in j may be summed up. Intro- 
ducing 

K(^) = j:Mn-,2j^-^^{-Zy (13) 
j>Q J-"'- 

X h{x 



ox (1 + x^ z 



\n+l ' 



we may then follow another time the reasoning of Section |lllj| replacing M„, Eq. (^), by 
M^{z) as given by Eq. ([T3|). We notice in passing that the modified moments M'^{z) are 
analytic functions of z throughout the complex plane except at z = —1. 

For simplicity we restrict ourselves to physical masses, i.e., to the positive real axis of z 
where the above integral representation is well defined. Then, two cases show up depending 
on whether rrih < 2eo or not. The former case is the one relevant to phenomenology but we 
consider both cases in turn for completeness. 

For TTifi < 2eo, we restore z = m\/{Ae^) and change the variable to x' = x/(l + 
x'^m1/{4el)). Then we have 



ML 



m 



2\ Al+ml/i^el)]-^ dx' h{x 



X 



"V4eg; Jo x' l-x^ml/{Aeiy 

with X understood as a function of x' . In this form we can easily follow the reasoning of 



Section [II B| because the integration range does not play a role until one cuts the amplitude. 
This cutting imposes x' = e^/X and thus results in a non vanishing imaginary part for 
A > Ao = Co + ml/ (4eo). Above this threshold 



ImM'iX) 



TT 3 2 /i(eo/A- 



— CtnC 



2^° ° l-m2/(4A2 



2"° ° JX 



/i(eo/A+), (14) 



where A+ = (^A + \JX'^ — rn?^/2. Dividing Eq. (|T4D by the fiux factor we obtain for the total 
cross section 

XI ri 



q-$fe(A) = 2 / dxG{x)a^g{xX+). (15) 

A^ — JO 
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Some comments are in order. We first stress that, as in the mh = case, one arrives 
at a simple partonic form of the cross section. This means, in particular, that the physical 
discussion we gave after Eq. (0), based on the subprocess ^ + g ^ Q + Q, still holds in 
the massive target case. Apart from the prefactor, the only modification between Eq. ( p!5[ ) 
and Eq. (|^ is the modification of the expression for the gluon energy in the partonic cross 
section o"$g. The change from xX to xX+ may be given a heuristic interpretation in the 
parton context using light cone coordinates^. In the $ rest frame let us choose the third 
axis along the hadron h momentum and form p+ = (A + \JX? — m\) / \/2 and p~ = m\/ (2p+) . 
In the parton picture the gluon causing the dissociation is picked up from the hadron h and 
has a negligible transverse momentum, and hence a negligible minus momentum. Its energy 
is then easily expressed in terms of x, the light cone (plus) momentum fraction of the gluon, 
and reads uj = xp'^/\/2, i.e., co = xX+. 

Our understanding of the prefactor is more formal. We observe that Eq. may be 
rewritten as 

ImM'iX) dX = lmM{X+) dX+, 
a relation which entails the (formal) identity 

MU{z) = M.i, Vz, (16) 

which can also be obtained from a direct comparison between Eqs. (0) and (§). 

Next, we point out that, as expected, the ruh 7^ corrections are sizeable only for small 
energies. The first aspect of these corrections is that, as above mentioned, the threshold is 
now located at 

Ao = eo + ^. (17) 

As in the massless case this corresponds to the need to find in /i a gluon with an energy 
UJ > eo sufficient to dissociate the $. With u = xA+ and x < 1 this gives a A-|_ threshold 
A+o = ^0; leading to Eq. (|1^. The second aspect is that the cross section behavior for 
A ^ Ao is given by Eq. (|TD|) with an argument A+ instead of A and a prefactor 

^+ / Co \ 



X^-ml \eo-ml/{4eo)J 

Let us now investigate the rrih > 2eo case. One may perform the same change of variable 
in the intervals [0,2eo/m/i] and [2eo/r?7,/j, 1] leading to 



MUmi/iAet,)) = —x 







x' 1 — ml/ (4eo) 



reo/mh dx' h{x) 



'[i+m2/(4.2)]-i x' x2m|/(4eg) - 1' 
The threshold becomes Aq = mh and one finds 

Im M'(X) = "^alel ( ^^^^^ + 4t^) • («) 



^ Such a connection between light cone variables and mh 7^ correction is discussed for deep inelastic 
scattering in ||l^ . 
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with \± = (-^ =t y^A^ — m|j /2. We notice that the first contribution is the one already 
obtained in the case < 2eQ. The second term is new but contributes only in the range 
< A < eo + m|/(4eo). We further point out that the relevant energy variable is now half 
the difference between energy and momenta, instead of half the sum for the first term, and 
that the present result Eq. (|18D is again consistent with the M_i sum rule ([T6|) . Finally, in 
the neighborhood of threshold the cross section now amounts to 



71 3 h{2€o/mh) 



2(AVm| 



1)2 



ao^o- 



(19) 



III. PHENOMENOLOGY 

A. Choice of parameters 

In view to give numerical values for leading twist total cross sections and thresholds, 
it is necessary to fix, on the one hand, the heavy quark mass mg and the quarkonium 
Rydberg energy eo, and, on the other hand, the gluon density in the target. 

1. Quarkonium sector 

The above described QCD analysis assumes that the QQ binding potential is well ap- 
proximated by the Coulomb part of the QCD potential ||^. Treating the IS and 2S heavy 
quarkonia as Coulombic states leads to 



Mqq{IS) 
Mqq{2S) 



2mQ - 
2mo - 



eo, 
fo 
A'' 



(20) 
(21) 



that is 



rriQ 



AMqq{2S)-Mqq{IS) 
Mqq{2S) - Mqq{IS) 



This gives for charm and bottom respectively (set (i)) 



eoc 
eofe 



0.78 GeV, 
0.75 GeV, 



rric 
rUb 



1.94 GeV, 
5.10 GeV. 



One way of estimating the applicability of the heavy quark analysis to charmonia and 
bottomonia is to compare the size of each Coulomb-state to typical confining distances. 
One may first evaluate the Bohr radius = 1/ y^eoniq, this gives aoc = 0.16 fm for charm 
and aob = 0.10 fm for bottom. Recalling that the IS-state root mean square is given by 
r{lS) = \^aQ one finds that the IS-state size remains somewhat below typical confining 
distances. We then consider that the LT analysis may be at least indicative of the behavior 
of IS-state cross section. Computing the size of 2S-states with r{2S) = ^/30ao, one sees that 
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the situation is much less favorable for 2S-states, especially for charmonium. The application 
of the framework to T' is given in Appendix 

In addition to the question of the validity of the computation of 2S-state cross section 
within the LT analysis this also led us to reconsider the above choice of parameters. For 



this we drop Eq. (21) and propose to fix the Rydberg energy to the energy-gap between 
the IS-state and the open flavor production. With Eq. ( pOf ) this is equivalent to putting 
"Trie = rrLD for charm and rUb = rriB for bottom. We then have the alternative set (set (ii)) 

eoc = 0.62 GeV, = 1.86 GeV, 
eob = 1.10 GeV, = 5.28 GeV. 

As we shall see, the cross section is only sizeable at large energy. In this region the magnitude 
of the cross section is driven by the factor OqCo = 1/ ^mgCo (see Eq ([TT|)). Taking set (ii) 
instead of set (i) results in a 20% cross section increase for charm and in a 22% cross section 
decrease for bottom. It turns out that this uncertainty is smaller than the one coming from 
scale fixing discussed in the next section. We will therefore limit our further considerations 
to cross sections obtained using set (i). 



2. Gluon distributions 

The other important input for the computation is the gluon density Gh in the hadron h 
considered. At this point one should remember that this density depends on the factorization 
scale /i (see paragraph following Eq. (P). In the <^-h cross section, not only Gh is a function 
of /X, but also a^g. Part of the dependence stands in the explicit coupling of the gluon with 
the $ constituents, corresponding to the factor as in a^e^ oc a^OQ, though the knowledge 
of the full dependence requires a complete one-loop calculation. Lacking such an analysis 
led us first to restrict ourselves to the so-called leading order (LO) gluon distributions (i.e., 
their evolutions are computed to one-loop) and correspondingly to the LO running coupling. 
Second, we investigated with some care the variation of our results with the choice of different 
factorization scales. 

We started this analysis with the prescription suggested in Ref. that ~ eo. For 
such a low scale the only available parameterization for the proton is that of Gliick-Reya- 



Vogt |16], labelled GRV94 LO^. These authors have proposed a parameterization for the 



pion too [ITj]. Concerning the x dependence of the various gluon densities, one should notice 
that the intermediate-a: region is rather well monitored, while the small-x region is poorly 
understood, especially for the pion case. Thanks to the HERA measurements the situation 
for the proton is much better. Including these data Gliick-Reya-Vogt have provided a new 
parameterization [|l^ but its lower scale is larger than that needed for this study. However 
we checked that at a scale large enough for both to be compared the difference between the 



gluon distribution of Ref. [|T8| and that of Ref. |]T6| is not significant 



We also examined, at a larger factorization scale, the consequences of different parame- 



terizations of the gluon distribution in the proton, considering in turn MRST98 LO ||T9[ and 
CTEQ5L g. 



^ It should be borne in mind that at such a low normalization scale the gluon distribution is less constrained 
by experimental data than by model assumptions. 
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B. Cross section variation with energy including target mass corrections 

As we shall see in the next section, the cross section at a given energy depends on the 
choice of parameters. Its general trend, however, is rather independent of a specific choice. 
We therefore begin our phenomenological study by discussing in this section those aspects 
that only weakly depend on the quarkonium parameters and gluon distributions. 

In Sect. |1| we have seen that one can distinguish two extreme energy regimes in the 
cross section: a threshold region and a high energy regime. The high-energy cross section is 
independent of the target mass and is given by 

o-(A)~a,,(A) = C(A/eo)' 

for gluon densities G{x) ~ const. /a; ^"'"'^ at small x. The constant C depends on details of the 
gluon density and on the parameters describing the charmonium sector (see, e.g., Eq. (pJ])). 
In every case discussed in Sect. [Ill A| the relative difference between this asymptotic cross 



section a^s and the full result is less than 25% for A > 30 eo. For J/ijj this translates roughly 
into ^/s>13 GeV and for T into > 23 GeV. 

For rrih = the cross section is very small in the threshold region. We found that it is 
less than 10% of for A < 2-3 eo, i.e., < 5 GeV for J/ip and 12 GeV for T. 

Let us now investigate the effect of a finite irih- Figure |l| shows the ratio of (T$/i(A,m/i) 
and (j$/i(A,0) for two hypothetical hadron masses such that m/j < 2eo, namely rrih = eo/5 
( "pion-like" , dotted) and ruh = eo. The latter case has been computed for a pion gluon 
distribution ("rho-like", dashed) as well as for a proton gluon distribution ("proton-like", 
solid). This ratio is plotted as a function of A/eo and is then identical for charm and bottom 
mesons. 

In addition to the shift of the threshold Eq. (P^, we observe that the inclusion of finite 
mass correction reduces the cross section close to threshold. This result is opposite to what 



is found in Ref. [|T0|^5 where it is argued that the target mass correction tends to increase 
the J/ip-p cross section near threshold. We notice that the mass correction is important 
only for A < 2-3 eo. This implies that it is of limited phenomenological interest since, as we 
have seen above, the cross section is very small in this low energy region. 

In Sect. II C we identified a different behavior in the case of heavy targets {rrih > 2eo). 



Figure |I] shows mass correction for a hadron with mass rrih = 3 eo and a gluon distribution Gp 
given by Ref. |l^ {dash- dotted). The cross section diverges at threshold (Ao/eo = rrih/eo = 
3), as can be seen in Eq. (0). The window for which the cross section gets sizeable is very 
narrow, however. Thus this threshold behavior has probably very little phenomenological 
implications. 



C. ^-h absolute cross sections 

1. Factorization scale dependence 

Before addressing in the next section the magnitude of the cross sections, the influence 
of the factorization scale is quantitatively investigated. More specifically the three choices: 
/i^ = eg, 2eQ, and 4 eg were considered. 



We give ample details on the comparison between the present approach and that of Ref. in Appendix ^. 
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FIG. 1: Ratio of the corrected ((j$/i(A, m/j)) over the uncorrected ((T$/i(A, 0)) cross sections as a 
function of A/eo- Calculations are performed for hadron masses nih = eo/5, eo and mh = eo, 3eo 
using the gluon density Gh{x,iJ. = 0.75 GeV) in the pion and in the proton respectively 
(see text). 



We made this investigation for the bottom channel with the parameter set (i). In addition 
to the heavy quark mass, this choice benefits from the numerical coincidence that® 




as expected for a large enough eo if the factorization scale is /i^ = Cq. 

Let us first define the T-p cross section, cryp~''°\ computed with the prescription /i^ = e^, 
to be that given by Eq. with the parameter set (i) for the bottomonium and the proton 

distribution GRV94 LO evaluated at /i^ = eg. We next define a^p^ at another scale to be 
that computed with GRV94 LO evaluated at fi^ and multiplied by the factor 

as{fJ,'^)/as{el), 

which takes into account the change of the coupling of the gluon to the $ constituents, and 
decreases down to 63% for /x^ = 4 eg. 

On Figure is shown the energy dependence of the T-p cross sections evaluated with the 
GRV94 LO gluon distribution at scales /x^ = (solid), 2 eg (dashed), and 4 eg (dotted). To 
bypass the question of the scale dependence of eo we restricted our study to the massless 
target case and studied the cross section as a function of A/eo. 



® For consistency we use here and in the foUowing the one- loop running coupling with = 3 and with the 
QCD scale determined by Refs. |l|l: A^^) = 232 MeV. 
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; : 




( 2 \ 

FIG. 2: Absolute cross section as a function of A/eo for different factorization scale: = 
(solid), ji^ = 2eQ {dashed), and fjp' = Ae^ {dotted). The gluon distribution used is GRV94 LO. 



We first remark that the higher the scale /i^, the larger (resp. smaller) the cross section 
at high (resp. low) incident energy. At high energy (A/eo ~ lO'^) the uncertainty may be 
as high as a 100%. The situation is much better in the range A/eo ~ 20 — 100 that is 
particularly relevant for phenomenology. We also notice that the running of as and that of 
Gp tends to somewhat compensate each other for the cross section at high energy. 

We also investigated the consequences of changing the parameterization of the gluon dis- 
tributions. Three leading order sets exist which may be evaluated at the scale /i^ = 4 eg. 
Figure displays the energy dependence of the absolute cross section a-^p using GRV94 LO 
{solid), MRST98 LO {dashed) and CTEQ5L {dotted) gluon densities. The energy depen- 
dence proves to be rather independent of a specific choice for energies A < 300 eo. At larger 
energies there is a rather strong dependence which leads to an uncertainty on the cross 
section comparable to that due to the scale variation (compare Figs. ^ and 0). The origin 
of this uncertainty is the poor knowledge of the gluon distribution at very low x. 



2. Cross sections using GRV gluon distributions 

We now turn to the discussion of the magnitude of the cross section. The J/ip and T 
cross sections are displayed on Fig. ^. They have been computed using the parameter set (i) 
with the GRV gluon density for the proton |16[ and the pion evaluated at /i^ = el 



These cross sections aj/^h i^^ft) and ayh {right) are found to strongly increase up to 
about 1 mb and 0.2 mb, respectively. The transition between low and high energy is situated 
around ^/sj/^h = 8 GeV and y/srh = 15 GeV, respectively. We also notice that, depending 
on the set of parameters chosen (resp. (i) and (ii)), the ratio Cj/^h/o'rh at high energy lies 
in the range 4-6, i.e., roughly the charm to bottom ratio of a^eo. 

The energy dependence of cr$7r turns out to be remarkably similar to the one in the 
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FIG. 3: Absolute crxp cross section as a function of A/eg using GRV94 LO (solid), MRST98 LO 
(dashed), and CTEQ5L (dotted) gluon distribution. The scale has been fixed for all distributions 
to /x^ = 4 en. 



proton channel, with a slightly smaller magnitude. This similarity is intimately related to 
the analogy that exists between the proton and the pion distributions in the GRV approach. 
Needless to say that lacking small x experiments for the pion it has not been possible to 
verify this analogy so far. With GRV distributions and at high energy (^s^h = 200 GeV) 
the ratio cr$,r/cr$p ~ 0.6, independent of the quarkonium considered. 

IV. CONCLUSION 

The operator product expansion analysis has been widely used in the analysis of deep 
inelastic scattering. Subsequently, these very techniques proved useful to investigate heavy 
quark systems p, |[ allowing the calculation of the J/ ip-p cross section within perturbative 
QCD. Such a cross section is of seminal importance in the context of heavy ion collisions. 

The present study is a continuation of the work of Bhanot and Peskin and of a more 
recent paper by Kharzeev and collaborators |jlO[|. Let us gather what have been carried out 



here. 

First, the leading twist forward scattering amplitude has been given a simple integral 
expression, entailing a partonic representation for the total cross section. Such a description 
had been found in a different way in Ref. 0] in the case of massless targets. 

Secondly, finite target mass corrections have been systematically incorporated. We 
showed that the cross section still assumes a partonic form though in terms of a modi- 
fied energy variable. In addition to a shift of the reaction threshold, finite mass corrections 
add to the suppression of the cross section at low relative energy but become insignificant 
far above threshold. In the case of heavy targets, however, we noticed that the cross section 
becomes large just above (and even diverges at) threshold. 
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FIG. 4: Absolute cross sections a^h as a function of the incident energy for J/ip (left) and T 
(right) with proton (solid), and vr (dashed). The gluon distributions Gp(x) and Gt^(x) used come 
from Ref. [16| and Ref. respectively, evaluated at iJ? = e^. 



Last, the energy dependence of aj/^h and arh has been investigated for several targets. 
We found that cr$,r and aq,p are strongly suppressed in the vicinity of the threshold. At 
large energy, the cross section is proportional to for a target with G(x) ~ const. /x^"^*^ at 
small X. With GRV gluon distributions this leads to slowly rising cross sections for both 
$-7r and However, we should emphasize that the small-x gluon distribution are not 

much constrained, especially that of the pion. Indeed the weak control we have on the gluon 
distribution, because of both the just mentioned poor small-x knowledge and the sizeable 
scale dependence, turned out to be the main source of uncertainty in the present approach. 
The $' cross section have also been investigated, although the relevance of a perturbative 
approach is not fully satisfied for 2S states. 

In addition to the hadron mass corrections considered here and beyond perturbative 
corrections, higher twist corrections may also lead to substantial modification for a not so 
heavy quark such as the charm. Indeed we noticed in Sect. [II B| that the threshold location 
may vary when one incorporates higher twist correction, making finite a cross section which 
is zero at leading twist. Considerations of this type of corrections is clearly outside the scope 
of the present study. Part of these corrections may be associated to the $ sector and in 
particular to the confining part of the heavy quark potential. Other corrections involving 
higher twist operators in the hadron target have presumably their counterparts in deep 
inelastic scattering. 

As compared to the other approaches mentioned in the introduction the J/ip-Tt cross 
section is very tiny at ^/s = 4 GeV. In the 4-7 GeV energy range it strongly increases driven 
by the intermediate x region in the gluon distribution. This x region is fairly under control 
thanks to the momentum sum rule and consequently the prediction, within the present 
approach, of cross section smaller than 1 mb below 7 GeV is rather robust. 

A consequence of such a small cross section at small relative energies is that destruction 
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of J/ip^s by comovers become very unlikely. We are presently studying whether our results 
allows already for an answer to the question whether a quark gluon plasma is formed in 
ultra-relativistic heavy ion reactions. Our present approach is, however, limited to the J/ijj 
cross section, since for the excited states the binding energy is not large as compared to the 
confining scale. 



APPENDIX A: SUM RULES FOR cr$ 



In this first appendix, we establish the sum rules satisfied by the leading twist total cross 
section aq,h- This allows us to make contact with a similar derivation done in the case 
rrifi = 00 and rrih 7^ [|lOl- The starting point is the expression for the leading twist 
amplitude (0). Expressing M„ = dnAn and recalling that An is the n-th moment of the 
gluon density in the hadron target (see Sect. |II A|) , we may write the series as 

M'{\m,) = alelj:iX/eor' f -x'^' G{x) 
k'>i " 

^{2k' + i)\^ ( ml ^\ 



Using the identities 

V ^ — \ n + 2i + -,-]z^ = B\n + 

3^-n\ \ 2'2) V 2' 2 

fb n 1 n 5 n n 

X 3-F2 - + -, - + -, 1 + n: - + -, 3 + -: z 
V4 24 2 2 2 2 

X 3F2 1, -, — : -, 4: z 
V 4 4 2 

the sums over j in Eq. (|A1|) leads to 

M'iX,mH) = alelJ2d2k'A{2k'){\/eof''' 

k'>l 

,2 



al€l(^jd2j dxxG{x) (A2) 



-0 

^ / 9 11 7 , m? 2' 
XsfJi, 4;-^.^^ 



where 



Ain) = C —x'^Gix 
Jo X 



^ f 5 n 7 n 5 n n ml ,^ 

X 3-F2 - + -, - + -, 1 + n: - + -, 3 + -: 1 x^ 

^^V4 2'4 2' '2 2' 2' 4eg ^ 
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We point out that Eq. (|A2|) is equivalent to the equation (15) of with shghtly different 
notations. 

As noted in performing the integral 

dX 



around a counterclockwise contour enclosing the origin gives the coefficient of A^^ in the 
amplitude Wrapping the contour around the u and s channel cuts, the contour at 

infinity giving no contribution^, 

2 r+oo 

- / d\ \'^^-Hm M'iX, rrih) = aleld2eA{2i) tf^. 

We put nih as a lower bound of the integral but this does not presume of the exact location 
of the threshold Aq (necessarily greater than or equal to rrih) implicitly contained in the set 
of sum rules (see below). Using the optical theorem one finally gets the sum rules for the 
cross section 

/ dX A-2^-VA2-m2a$,(A) = -alel d2iA{2i) ef', (A3) 

which is what Bhanot and Peskin found (Eq. (3.10) in Ref. 0]) in the limit of massless 
target. In order to compare these results with |T^, we introduce the variable y = rrih/X to 
get 

''dyy^'~^^l^^a^^{m^/y) = {'^V'' I{2^)A{2^), (A4) 



with 



I{n) = -aleo d^, 



as introduced in [llOl. The relation (|A4|) gives the set of sum rules that should replace that 



given by Eq. (16) of Ref. ||T0[, where the prefactor {rrth/ eof'^'^ is missing. 

In order to verify that the sum-rule formalism leads to the same result as we obtained 
in the present study we need to solve Eq. ( |A4|) for o"$/j. This may be done using Laplace 
transform techniques. Eq. ( |A4| ) writes 

dyy'^'^yi - y^(T^h{mh/y) = g{n), n = 2,4, . . . , 

^ 
with 

n-l 



gin) = (^)"' I{n)A{n). 



eo / 

Defining x = —Iny and n' = n — 1, we get the relation 



r+co I 

/ dx exp(— n'x) \/l — exp(— 2x)(T$/i(m/i exp(x)) 
JO v! / 



=/(^) 

g{n' + l), n' = l,3,.... 



^ Both the singularity pattern and asymptotic behavior of the scattering amphtude, which are necessary 
for the present construction to hold, are best studied within the approach of Sect. ||. In Ref. [Q, these 
properties where assumed from general properties of the elastic scattering amplitude. 
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that is g{n' + 1) is the Laplace transform of f{x). 

The above relation can be uniquely extended to a complex argument u (at least for 
Rez/ > 6, see below) instead of the listed n'. We can then obtain f{x) by inverting the 
Laplace transform 

fix) = 77— / duexp{ux)g{u + 1), (A5) 

where u'^ is an arbitrary real chosen so that the integration contour is located "at the right" of 
all singularities of g{i' + 1). For gluon distribution behaving as x~^^^^'' for x — 0, integrals 
of the type Jq dx x''G{x)3F2 are well behaved provided Reu > 6. As a consequence we 
performed the inverse Laplace transform choosing u'q > 6. 

Setting u = Uq + iu = uq ~ 1 + iu, and then A = rrih exp(x), we thus obtain 

yr^^exp(^^2x) 0-$ (m/, exp (x) ) 

1 f+oo 

du exp[{uo — 1 + iu)x]g{i'o + iu), 



27T 



and then 



o-^hW = , / du 



-00 




27r^A2 - m 



2 J— 00 

h 




I{uo + iu)A{iyo + iu), (A6) 

for every Uq > 1 + 6. Let us remark that the eo energy scale appearing in the first factor 
of the integrand would be rrih with the sum rules proposed in |TD[. The difference between 



these two results is important when one considers the m/j — > limit, in which case the latter 
expression is ill-defined contrarily to ours, given in Eq. (|A6|) . 

We carried out a numerical evaluation and found out that it reproduces the results ob- 
tained in the main body of the paper. One critical point in the comparison between the two 
approaches is the verification that the threshold is located at the predicted value, that is, in 
particular, that the numerical result is compatible with zero below this threshold (a point 
which is far from evident when one looks at Eqs. or (|A3|) ) . 



APPENDIX B: 2S STATES 

In this appendix, calculations of the $ cross sections are extended to the 2S states 
The modification amounts to replacing the IS coefficients d^^'^^ (denoted dn so far) by the 
2S coefficients 

16^ r(ra + 5/2)r(5/2), . 

'S'-sn! r(„ + 7) <^°" + 
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Expressing d^^^ function of the IS coefficients 

allows one to get an integral representation for the 2S coefficients 



42^) =4" [' ^xy'^''\x), (Bl) 

JO X 



with 

f^''\x) = ^^x'/\l-x)n^-Qxr, 



c 



which is the ingredient needed to carry out the procedure outlined in Sect. ||. The changes 
are that the function h has to be replaced by h^"^^^ = G f^"^^^ and that it is now evaluated 
at e/A,e/A+, or e/A_ instead of eo/A,eo/A+, or eo/A_ where e = eo/4 is the binding energy 
of the 2S state. The partonic expression is thus similar with 

= 16^aoeo Oiuj - e). 

The energy dependence cr^'pis) has been computed using the parameter set (i). The cross 
section diverges at threshold (^i 11 GeV) and decreases to a minimum of about 4 mb at 
a center of mass energy 0.6 GeV above threshold. Then the cross section increases smoothly 
((Tt'p(s) oc s^) and reaches 30 mb by y/s ~ 200 GeV. We notice that at high energy the ratio 
(T$//i/(T$/i ~ 20 for all incident hadrons. Since r{2S) = vTOr(lS'), this result lies somewhat 
above the geometrical expectation r^(2S')/r^(15'). 

We have already insisted on the fact that the LT perturbative analysis is most likely not 
adequate to describe the ip' channel. In this case, the cross section amounts to o^q,rp{^/s) = 
30 mb at = 10 GeV. 
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